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Abstract. A general theory of homeoidally striated density profiles where no di- 
vergence occurs, is adapted to cuspy density profiles, with a suitable choice of 
, the scaling density and the scaling radius. A general formulation of some physi- 

o . 

calculated in connection with a special class of rotational velocity fields), inertia 

tensor, and self potential-energy tensor, is performed. Other potential-energy ten- 

sors involving two density profiles where the boundaries are similar and similarly 

q . placed, are also expressed. Explicit results are attained for three special cases of 

5— i ' 

. physical interest: NFW (e.g., Navarro et al. 1997) and MOA (e.g., Moore et al. 

CZJ . 

1999) density profiles, which fit to a good extent the results of high-resolution sim- 
ulations for dark matter haloes, and H (Hernquist 1990) density profiles, which 
closely approximate the de Vaucouleurs r 1 / 4 law for elliptical galaxies. The virial 

;_h ' 

theorem in tensor form for two-component systems is written for each subsystem, 
and applied to giant elliptical galaxies. The predicted velocity dispersion along 
the line of sight, in the limiting case where a principal axis points towards the 
observer, is found to be consistent with observations except for (intrinsic) E7 
configurations where the major axis points towards the observer. If dark matter 
haloes host an amount of undetected baryons about twice as massive as the stellar 
subsystem, and undetected baryons trace non baryonic matter therein, two main 
consequences arise, namely (i) velocity dispersions along the line of sight are lower 
than in absence of undetected baryons, and (ii) dark matter haloes are dynam- 
ically "hotter" than stellar ellipsoids, the transition occurring when the amount 
of undetected baryons is about one and a half times that of the stellar subsys- 
tem. In this view, both the observation that the temperature of the extended hot 
gas exceeds the central stellar temperature, and the fact that the non baryonic 
matter is dynamically "hotter" than the stars, are a reflection of the presence of 
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undetected baryons, which trace the dark halo and are about twice as massive as 
the stellar ellipsoid. 



Key words. Cosmology: dark matter - galaxies: haloes. 
1. Introduction 

According to standard CDM or ACDM cosmological scenarios, large-scale celestial ob- 
jects such as galaxies and clusters of galaxies, are made of at least two components: one, 
baryonic and more concentrated, embedded within one other, non baryonic and dissipa- 
tionless, usually named dark matter halo. After a wide number of both analytical and 
numerical studies (e.g., Cole & Lacey 1996; Navarro et al. 1995, 1996, 1997; Moore et al. 
1998, 1999; Fukushige & Makino 2001; Klypin et al. 2001), it has been realized that dark 
matter haloes which virialize from hierarchical clustering show universal density profiles, 
p = p(r; p\r'), where p^ is a scaling density and r' is a scaling radius. In this view, 
smaller haloes formed first from initial density perturbations and then merged with each 
other, or were tidally disrupted from previously formed mergers, to become larger haloes. 

The density profile is (i) self-similar, in the sense that it has the same expression, 
independent of time (e.g., Fukushige & Makino 2001), and (ii) universal, in the sense 
that it has the same expression, independent of halo mass, initial density perturbation 
spectrum, or value of cosmological parameters (e.g., Navarro et al. 1997; Fukushige & 
Makino 2001). A satisfactory fit to the results of numerical simulations is the family of 
density profiles (e.g., Hernquist 1990; Zhao 1996): 



for a suitable choice of exponents, a, (3, and 7. 

This family includes both cuspy profiles first proposed by Navarro et al. (1995, 1996, 
1997), (a,f3, 7) = (1,3,1), and the so called modified isothermal profile, (a, /?, 7) = 
(2,2,0), which is the most widely used model for the halo density distribution in anal- 
yses of observed rotation curves. It also includes the perfect ellipsoid (e.g., de Zeeuw 
1985), (a, P, 7) = (2,4,0), which is the sole (known) ellipsoidal density profile where a 
test particle admits three global integrals of motion. Finally, it includes the Hernquist 
(1990) density profile, (a, /?, 7) = (1, 4, 1), which closely approximates the de Vaucouleurs 
r 1 / 4 law for elliptical galaxies. In dealing with the formation of dark matter haloes from 
hierarchical clustering in both CDM and ACDM scenarios, recent high-resolution simu- 
lations allow (a,/3, 7) — (3/2,3,3/2), as a best fit (e.g., Ghigna et al. 2000; Fukushige & 
Makino 2001; Klypin et al. 2001), as first advocated by Moore et al. (1998, 1999)0. 

Send offprint requests to: R. Caimmi 
1 More precisely, an exponent a — 1.4 was derived by Moore et al. (1998), while the value 

a = 1.5 was established by Moore et al. (1999). 
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Though Eq. ([!]) implies null density at infinite radius, the mass distribution has nec- 
essarily to be truncated for two types of reasons. First, the presence of neighbouring 
systems makes the tidal radius an upper limit. On the other hand, isolated, over-dense 
objects cannot extend outside the Hubble sphere of equal mass. Second, the total mass, 
deduced from Eq. (|l]) for an infinitely extended configuration, is divergent, at least with 
regard to the special choices of exponents, (a, f3, 7) = (1, 3, 1), hereafter quoted as NFW 
density profile, and {ot,j3, 7) = (3/2,3,3/2), hereafter quoted as MOA density profile. 
The region enclosed within the truncation boundary has to be intended as representa- 
tive of the quasi static halo interior, leaving aside the surrounding material which is 
still infalling. It is worth remembering that the total mass can be finite even for an in- 
finitely extended configuration, provided the related density profile is sufficiently steep 
e.g., (a, (3, 7) = (1, 4, 1), hereafter quoted as H density profile. 

In dealing with numerical simulations of dark matter haloes, it is usual to take into 
consideration spherically averaged density profiles (e.g., Cole & Lacey 1996; Navarro et al. 
1997; Fukushigc & Makino 2001; Klypin ct al. 2001), or in other terms spherical isopyenic 
(i.e. of equal density) surfaces. On the other hand, spin growth by tidal interactions 
with neighbouring objects, in expanding density perturbations, demands ellipsoidally 
averaged density profiles (e.g., Doroshkevic 1970; White 1984; Mailer et al. 2002; Jing 
& Suto 2002), or in other terms ellipsoidal isopyenic surfaces. As a best compromise 
between intrinsic simplicity and unavoidable necessity, our attention will be devoted to 
homeoidally striated configurations, i.e. the isopyenic surfaces are similar and similarly 
placed ellipsoids. 

Galaxies and cluster of galaxies may safely be idealized as two subsystems which link 
only via gravitational interaction, in such a way that each component is distorted by the 
tidal potential induced by the other. Then the application of the virial theorem in tensor 
form may be performed either to the whole system or to each subsystem separately. 
Towards this aim, the explicit expression of the potential-energy tensors are needed. 
Though some results are available in literature (e.g., Brosche et al. 1983; Caimmi & Secco 
1992; Caimmi 1993, 1995), the related density profiles exhibit no central divergence, or 
"cusp", in contradiction with Eq. ([!]) when 7 > 0. 

The present attempt aims mainly to (i) formulate a general theory of homeoidally 
striated density profiles with a central cusp; (ii) devote further investigation to a few 
special cases which are consistent with the results of both observations and simulations 
related to galaxies and cluster of galaxies; (iii) apply to galaxies, as the analogon to 
clusters of galaxies performed in an earlier paper (Caimmi, 2002). 

The current paper is organized in the following way. The general theory of homeoidally 
striated, density profiles with a central cusp, is performed in Sect.^J. The general results 
are particularized to NFW, MOA, and H density profiles, which provide good fits to the 
results from both observations (e.g., Geller et al. 1999; Rines et al. 2001) and numerical 
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simulations (e.g., Klypin et al. 2001), in Scct.||. An application to elliptical galaxies is 
performed in Sect.^j. Finally, some concluding remarks are reported in Sect.|^. Further 
details on fitting simulated and theoretical, self-similar, universal density profiles, are 
illustrated in the Appendix. 

2. Homeoidally striated density profiles with a central cusp 

A general theory of homeoidally striated density profiles, where no divergence in the 
density occurs, has been developed in an earlier approach (Caimmi 1993, hereafter quoted 
as C93). Here we take into consideration a more general approach, where a central cusp 
may occur. The isopyenic surfaces are defined by the following law: 

p = p t f(0 ; /(i) = i ; (2a) 

l=\ \ a l> 

where (ft = p(l) , al, are the density and the semiaxes, respectively, of a reference isopy- 
enic surface, and 3 corresponds to the truncation isopyenic surface, related to semiaxes 
ae- 

Then the scaling density, , and the scaling radius (related to the reference isopyenic 
surface), r\ correspond to a single boundary, which allows the description of cuspy 
density profiles. On the other hand, centrally regular density profiles defined in C93 fail 
in this respect, as the scaling density coincides with the central density and the scaling 
radius attains the external surface. 

The assumption that the system is homeoidally striated (Roberts 1962) implies the 
following properties (e.g., C93): 

a' da' , 

£=-f 5 d£=^ ; p= 1,2,3 ; (3a) 

dp dp 

; p = w ; S^ = ^ ; (3b) 

where a' are the semiaxes of the isopyenic surface under consideration. 



The generic point on an isopyenic surface obeys the equation: 

^ = 1; (4) 



r 2 sin 2 9 cos 2 6 r 2 sin 2 9 sin 2 i 



Uj-y LA3 Ujy u*2 1*3 

where P(ii, £2, £3) = P(r, 9, <f>). Taking spherical coordinates into consideration, Eq. ( 
may be cast into the equivalent form: 

sin 2 9 cos 2 cb + ^ s^ 2 9 sin 2 + ^ ^ cos 2 * 

(rt) 2 ( a t)2 a' 2 V ( rt ) 2 (4) 2 a ? 

( r t)2 f r U2 

in 2 9 cos 2 6 + -^-r— sin 2 9 sin 2 6 + ±-r— cos 2 1 



r 2 (rt 


2 (4) 2 


^) 2 (a\) 2 a? 


r 2 1 


[(rt) 2 ^ 


(rt)2 £2 


(4)2 ^ 



1 ; (5) 

(4) 2 ' (4)° 

where r* is the radial coordinate of the point, P'(r',9, (f>), on the reference isopyenic 
surface. Accordingly, the expression within brackets in Eq. (j^) makes the particularization 
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of Eq. (Q) to the reference isopycnic surface, and necessarily equals unity. Then Eq. (||) 
reduces to: 

C = 4 ; (6) 

r I 



regardless of the radial direction. The particularization of Eq. (|6J) to the truncation bound- 
ary yields: 

S = 4; (7) 

regardless of the radial direction. 

The above results may be summarized as follows. Given a homeoidally striated den- 
sity profile, and two generic points with coinciding angular coordinates, P(r, 6, <j)) & n d 
P'(r\ 9, </>), placed on two isopycnic surfaces, r(9, 4>) and r'(9, </>), the ratio of the radial 
coordinates, £ = r/r\ does not depend on the radial direction. Then Eqs. ( pa] ) follow 
from (JsJ) as special cases. 

2.1. Mass and inertia tensor 

The volume, the mass, and the inertia tensor of an infinitely thin homeoid bounded by 
isopycnic surfaces, may be expressed as (e.g., C93): 

dS" = 47ra|44C 2 d£ = 35^ 2 d£ ; (8) 
AM' = 4V444/(0£ 2 d£ = 3A/t/(£)£ 2 d£ ; (9) 

A-jr 

dl' pq = — p ^ pq (alf a y r f(0C 4 dC = 5 pq MHal) 2 fm 4 d£ ; (10) 

where and are the volume and mass, respectively, of a homogeneous ellipsoid, 
bounded by the reference isopycnic surface, and with the same density, p', as at the 
reference isopycnic surface; and S pq is the Kronecker symbol. 
Let us define, according to Roberts (1962): 



= 2 / /(Of dr ; (ii) 



f 



from which the following relations are easily derived: 

F(E) = ; ^ = -2£/(0 ; (12) 



an integration by parts of Eq. ([LI]) shows that: 

f(OC^=^- T F(OC' 2 ^ ; n>l ; (13) 



o 2 Jo 

which allows the calculation of the total mass as: 

M = v mas M X ; (14a) 
3 



in fj.s 



2 



F(£)d£ ; (14b) 



(i 



M+ = ypMot t ; (14c) 
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and the inertia tensor as: 



Ipq = Sp q Ui nr M' (a' p ) 2 ; (15a) 

Vinr=\j~ F{Oe^ ; (15b) 

where the coefficients, v mas and Vinr, are shape-independent and may be conceived as 
profile factors (C93). 

The mass enclosed within a selected isopyenic surface, p = p'/(0) < £ < 2, is 
obtained by integration of Eq. (||) . The result is: 

M(£) = 3Mt f(Oe 2 de ; (16) 



^) = f^ = 3p^ ^ /(n£' 2 d£' ; (17) 



JO 

and the related mean density is: 

where Eqs. (||) and (|l6|) have been taken into consideration. 
2.2. Tensor potential and self-energy tensor 

The tensor potential induced by an infinitely thin homeoid bounded by isopyenic surfaces, 
on a generic point placed at its interior, is (e.g., C93): 

dV#**> = 27rG6 pqP \al) 2 A p f(^d( ; (18) 

where G is the constant of gravitation, p = 1, 2, 3, q = 1, 2, 3, and A p are shape factors 
which, ipso facto, depend on the axis ratios only (e.g., Caimmi 1992; therein defined as 
a p ). 

The tensor potential induced by all the homeoids enclosing a given point, on that 
point, is (e.g., C93): 

V&**> =^G5 mP \alfA p F{Q ; (19) 

where £ p denotes the isopyenic surface passing through the point, and Eq. ( pi] ) has been 
used. 

The interaction-energy tensor related to an infinitely thin homeoid bounded by isopy- 
enic surfaces, embedded into another infinitely thin homeoid of the same kind, is (C93): 

d 2 (E mt ) pq = ~^ 2 G5 pq {p^) 2 {alf a\alA p f{0if{Oi' 2 d^d^ 1 . (20) 

if the two boundaries coincide, and d 2 (E int ) pq (£=£')= d 2 (E se i) pq , i.e. the 

interaction-energy tensor reduces to the self-energy tensor. Bearing in mind that Eq. ( p0| ) 
remains unchanged by replacing £' with £ and vice versa, the self-energy tensor of the 
whole mass distribution may be written as (E se i) pq = 2 J J s d (Ei n t) pq ; the result is 
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(e.g., C93): 




P I 



(21b) 



(21a) 



(21c) 



where v ae \ is a profile factor, i.e. shape- independent, and e pq are axis ratios. 

Let P(r, #,</>) be a generic point on an isopyenic surface, p = p^f{£,), < £ < 5. A 
test particle at P is centrifugally supported with respect to the centre of mass, provided 
the following relation holds: 



where rj represents the deviation from the spherical limit and/or the Roche-like (i.e. 
a mass point surrounded by a masslcss, ellipsoidal atmosphere) limit, and v eq is the 
rotational velocity with respect to the centre of mass. Due to Newton's theorem, a generic 
homeoid bounded by isopyenic surfaces, related to £' > £, £" > £, exerts no action on P 
(e.g., Chandrasekhar 1969, Chap. 3, §17). In this respect, only the mass, M(£), enclosed 
within the isopyenic surface under consideration, is effective. 

The particularization of Eq. ( p2|) to a point, P^(r^, 6, (j)), on the reference isopyenic 
surface, reads: 



(rt)2 r t 
and the combination of Eqs. @, (^2|) and (|23|) yields: 

= 1 Af (Q y(r, 6, 0) 
vUri,0,4>) f M(l)»7(rt,fl,^) ' 



provided the point on the generic isopyenic surface, P(r,9,(p), and the point on the 
reference isopyenic surface, (r^ , , , are aligned along the same radial direction. 

It is worth noting that the density profiles of interest, i.e. decreasing from the centre 
to the boundary, lie between two extreme situations: the homogeneous limit and the 
Roche-like limit, where the ratio, r](r, 9, (j>)/rj(r' , 9, (f>), reduces to unity in both cases. 

If the angular coordinates are specified, Eq. ( pif ) takes the simpler form: 



GM(€)T,(r,0,<l>) 



vl q {r,9,cj>) 



(22) 



r 



<(£) _ 1 M(£) ?7(0 



(25a) 



{vl q f £M(1) ryt ' 
v l q = v eq{r\9,(j)) ; = r]{r\e,4>) 



(25b) 



where it is intended that the velocity profile corresponds to the (r\ 9, </>) direction. 
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2.3. Angular-momentum vector and rotational-energy tensor 

In dealing with angular momentum and rotational energy, the knowledge of the rotational 
velocity field is needed. Let us define the angular-momentum vector and the rotational- 
energy tensor, respectively, as: 

J, = j j J n p d 3 i u + J J J n p d 3 i qq 

Q e d 3 I pi + f f f Q q d 3 I pq ; p £ ; q ? £ ; (26) 
s J J Js 

2(E ro t)pq = J J J Q P Q q d 3 Iu — J J J QtQ p d 3 I q t 

i3 t i / / / r>2 j3 



n q n e d 3 I ep + / / / Qjd 6 I pq ; p^i ; q^£ ; (27) 
s J J Js 

where fi = (Qi,fi,2, ^3) is the angular-velocity vector, and d 3 I pq is the inertia tensor, 
related to a generic, infinitesimal mass element. 



The symmetry of the tensors under consideration, due to Eq. (15a), together with the 
conditions, p ^ £, q ^ £, implies the following values for the index, £: 



1 = 1 , pV q = 2 , p^lkq^l 
£=2 , p vg=3, p^2hq+2 
£ = 3 , pVq= 1 , p^3 & q^3 



(28a) 
(28b) 
(28c) 



where V and & are the semantic symbol of disjunction, "or" , and conjunction, "and" , 
respectively. 

The preservation of ellipsoidal shape imposes severe constraints on the rotational ve- 
locity field. Leaving an exhaustive investigation to more refined approaches, our attention 
will be limited to a restricted number of special cases, namely: (i) rigid rotation about a 
principal axis, and (ii) differential rotation about a symmetry axis, in connection with a 
particular class of velocity distributions. By "rotation" it is intended, of course, circular 
rotation. Accordingly, f2 = (<5i s f2i, S2 S ^2, <^3s^3) provided x s is assumed as rotation axis. 

The angular momentum and the rotational energy, related to axis x s , are: 

J(s) = J s = J J J^Q s d 3 I pp + J J J^ s d 3 I qq ■ p^q^s ; (29) 

Erot(s) — {E ro t)s = (E ro t)pp + (E ro t)qq 

n 2 s d 3 i pp + ([[ n 2 s d 3 i qq ■ p^q^s ■ (30) 

is J J Js 

where s — 3 is usually assumed, and then p = 1, q = 2, or vice versa. For the sake of 
brevity, from this point on J s and (E rot ) s shall be denoted as J and E rotl respectively. 



In the special case of rigid rotation about a principal axis, Eq. (|27|), by use of (15a), 
reads: 

{E r ot)pq = T^Ipq^l = ^ pq (l - S pr )v inr M f (a* ) 2 f^ = ^S pq (l - S pr ) v mr A/ 1 * {v\ ot ) 2 ; (31) 



R. Caimmi and C. Marmo: density profiles with a central cusp 9 
and Eqs. @, @, taking s = 3, read: 

J=(/ii+/ 22 )^ = ^r l rM t a t 1 (l + e2 1 )^ ot ; (32) 

S ro4 = i i / inr Mt(l + e| 1 )(4 ot ) 2 ; (33) 

where v rot — x r and wJot = (iVoi)i is the rotational velocity at the end of the major 
semiaxis, a\, with regard to the reference isopycnic surface. 

In the special case of differential rotation about a symmetry axis, our attention shall 
be restricted to (rotational) velocity distributions which satisfy the counterpart of Eq. (g), 
namely: 

Vrotjr, 9) _ ty ot (a',0) _ , . 

Vrot{r 1 ,0) v rot (a\0) 

where (r, ff), {r\9), represent a point on a generic and reference isopycnic surface, re- 
spectively, along the same radial direction, and (a',0), (cr,0), represent the end of the 



corresponding, major semiaxes. Due to Eqs. (3a) and (|6|), Eq. (34) may be written under 
the equivalent form: 

"M) = o) (v ~ 
fi(rt,0) n( t,o) ' 1 1 

where rotational velocity is replaced by angular velocity. It is worth noting that, in 

particular, either rotational or angular velocity is allowed to be constant everywhere. 

Both rotational and angular velocity are independent of the longitudinal angle, 4>, 

according to Eqs. (pi} ) and (|35|), respectively. It follows that sections of infinitely thin, 

homogeneous homeoids, by infinitely thin layers normal to the symmetry axis, do rotate 

rigidly. The related volume, inertia tensor, angular-momentum vector, rotational-energy 

tensor, are (e.g., Caimmi & Secco 1993): 

a 2^V)3^ si nfldedfl 

(cos 2 + sin 2 0)3/ 2 ' 1 ' 



d% = d% = ^ lg Mt(at)V/(^ 2 f?2 ded L 5/2 5 ( 37a ) 
4 (cos 2 9 + e 2 sin 9yl l 

,a T 3 P „, + / +N 9 4r/ , X ,4 sin 9 cos 2 9 df d9 ,„„ , 

d 2 J 3g = -6 3g MV V/ ^ ; (37b) 

2 (cos 2 # + e 2 sin t^/ 2 

d 2 J = f!(£,0) [d 2 /n+ d 2 / 22 ] ; (38) 

d 2 (E rot ) pq = ^S pq (l-S p3 )n 2 (C,9)d 2 I pq ; (39) 

where, due to axial symmetry, a\ = a\ = a\ e 2 i = 1, e 31 = e, and 9 is the azimuthal angle. 



The angular-momentum vector and the rotational-energy tensor, by use of Eqs. (14c) and 
(H)~(H)> take the equivalent form: 

^,)„ - |.„(i-^)M. (a *)^v(i)i'^«.o) g^j (m 4°;^ in d ; )5/2 ;(«) 
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where, due to Eqs. (||), @, the ratios, v rot {l,0)/v rot (l,O), Q(1,0)/Q(1,O), only depend 
on the azimuthal angle. 

An integration over the whole, homogeneous, infinitely thin homeoid yields: 



dJ = 2r? amn AfV) 2 O(£,0)/(£)£ 4 d£ ; (42) 

_ 3 4 f + ^ 2 n(l,g) sin 3 fldfl 
Vanm 4 £ J_ v/2 0(1,0) (cos 2 + e 2 sin 2 0)5/2 ' ^ 

d(£?rot)p<? = <^pg 

(l-« P s)i|™iW t (a , ) , fl J K 1 0)/K){ 4 d( ; (44) 
3 4 /" +7r/2 fi 2 (l,0) sin 3 0d(9 



??ro ' 8 £ i_ V2 fi 2 (l!o)(cos 2 + e 2 sin 2 0)5/ 2 ; (45) 
finally, an integration over the whole system yields, for the angular-momentum vector: 

J = 2rj anm i' anm M i< vl ot a { ; (46a) 

*— = j'm^) m(,Ai • (46b) 

and for the rotational-energy tensor: 

(E ro t)pq = S pq (l - 8p?,)ri rot v ro tM ] {vl ot f ; (47a) 

(47b » 

it is apparent that rj anm , rj rot , depend on the integration over the boundary and may be 
thought of as shape factors, even if they may also be shape- independent, whereas v anm , 
v r ot, are genuine profile factors. 

The special case of rigid rotation, implies the relation: 

the special case of constant rotational velocity on the equatorial plane, implies the rela- 
tion: 

lH = a l = I ■ (49) 
0(1, 0) a' £ ' [ ! 

the special case of constant rotational velocity everywhere, implies the relation: 

n(g,fl) _ a' _ (cos 2 6> + e 2 sin 2 fl) 1 / 2 
O(£,0) ~ rsin(9 ~ esintf ' 

the corresponding values of the shape parameters, rj anm , r] ro t, and profile parameters, 

Vanm, v rot , are listed in Table [l| 

2.4. Other potential-energy tensors for subsystems with similar boundaries 

The above results are related to a single subsystem: strictly speaking, all the quantities 
discussed in the last Section should be labelled by the index, u, in connection with uth 
subsystem (e.g., Caimmi & Secco 1992), but it has been omitted for the sake of clarity. 
On the other hand, the formulation of the interaction-energy tensor and the tidal-energy 



(50) 
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case RR RC CV 

isonvcnic surface - 1 "(€,9) _ i <w (£,<>) _ i 

isopycnic surtace £=2(sToJ - 1 ?T(c7)J ~~ 1 *v ot (€,o) - 1 

equatorial plane ^^^-^j - l „ rof (i ;0) - 1 „ roi (i,(]) - 1 

T]anm 1 1 37r/8 

J?rot 1/2 1/2 3/4 



atrm 



f« meat j\ f<LK''K 



Table 1. Values of the shape parameters, i] anm , r q rotl and profile parameters, v a n mi v rotl 
related to (i) homeoidally striated ellipsoids in rigid rotation about a principal axis (RR) , 
and (ii) homeoidally striated spheroids in differential rotation about the symmetry axis 
where the isopycnic surfaces are: either in rigid rotation and the rotational velocity on 
the equatorial plane is constant (RC), or in cylindric rotation and the rotational velocity 
is constant everywhere(CV). 



tensor necessarily involves at least two components. The related calculations turn out to 
be very difficult in the general case, and for this reason our attention shall be limited 
to subsystems with similar boundaries, and homeoidally striated density profiles with a 
central cusp. 

Using a similar procedure as in C93, where no (central) cusp occurs in density pro- 
files, but bearing in mind that the scaled density and the scaled radius have a different 
definition, yields the expression of the interaction-energy tensors: 
G(M^) 2 

[(Eij)int}pq = —Spq — ~ +~ {^i^int^^Ap J (51a) 

Mint = ~rrf [tu^Hi/yt) +™ (ext) (S i /2/ t )] ; (51b) 

int]pq — \{^ij")int\pq i 

(51c) 

the expression of the tidal-energy tensors: 
G(M^) 2 

{{Eij)tid]pq — —5 P q , ' — {vij)tid£ P 2tp3A p ; (52a) 

(«y)tw = -|»»t «;<"*) (3,/yt) ; (52b) 
G{m}) 2 

[(Eji)tid]pq = Spq , {Vji)tid£p2tp3A p ] (52c) 

(v ji )m = -l^ W {int \z i /y f ) ; (52d) 

and the expression of the residual-energy tensors: 
G(M^) 2 

— {vij) res e P 2e p zA p ; (53a) 

K T )i 

{^)res = ~m! [w^HEt/y^) - «;<*»*> (Hi/yt)] ; (53b) 

[{Eji^) res ]pq = [(Eij^) res ]pq J (53c) 



12 R. Caimmi and C. Marmo: density profiles with a central cusp 

where the functions, u/ mt ) and w^ ext \ are defined as: 
fV a p. 

W (™t) {v)= / . (54a) 

JO C K] 

rv A p. 

W (^( V )= / i&d^ ; (54b) 

Jo Q si 

the parameters, m) and y', are defined as: 

mt = -4 ; y f = ^ ; (55) 

M} r\ 
and the combination of Eqs. (J?]) and (|5^) yields: 

U = ^ ; y = f ; (56) 

finally, the virial-energy tensors (e.g., Caimmi & Secco 2002), [{E uv ) vir ] pq = [(-E , u ) se ;]p, 2 + 
[(E uv )ud]pq, take the expression: 

G(Mt) 2 

\{,E UV ^)vir\pq — Spq - r~ {yuv)vir^p1^-pZ^-p 5 (57a) 

(a;)i 

(Wut>)mr = (w«)aei + (^««)«d ! U = i,j \ V = j,i ] (57b) 



in conclusion, Eqs. (|5lj)-(57) allow the calculation of potential-energy tensors for 
homcoidally striated density profiles with a central cusp, related to two subsystems with 
similar and similarly placed boundaries. 

3. Special cases 

The physical parameters related to the density profiles expressed by Eq. (Q), and inves- 
tigated in Sect. 2, may be explicitly calculated, provided the exponents, {a, (3, 7), are 
specified. Our attention shall be limited to two special cases, which are consistent with 
high-resolution simulations of cold dark matter haloes (e.g., Klypin et al. 2001): NFW, 
first recognized by Navarro et al. (1995, 1996, 1997), and MOA, first recognized by Moore 
et al. (1998, 1999). As far as we know, no method has been clearly outlined for fitting 
simulated and theoretical, self-similar, universal density profiles. A possible procedure is 
sketched out in Appendix^ for a single matter distribution which obeys Eq. (|l|), and in 
Appendix^ for two different mass distributions which obey Eq. (Q). In addition, a third 
density profile will be taken into consideration, which is consistent (together with NFW) 
with observations (e.g., Geller et al. 1999; Rines et al. 2001; Holley-Bockelmann et al. 
2001): H, first proposed by Hernquist (1990). 

With regard to the shape parameters, r) anm , rj rotl and the profile parameters, v anm , 
v ro t , it is worth remembering that they are listed in Table |] for the velocity profiles dis- 
cussed in Sect.^j, where differential rotation corresponds to axisymmetric configurations 
with constant rotational velocity, either on the equatorial plane or everywhere. 



R. Caimmi and C. Marmo: density profiles with a central cusp 
3.1. The NFW density profile 



13 



The NFW density profile corresponds to the choice of exponents (a,/3, 7) = (1,3,1) in 
Eq. ([|) . Accordingly, Eq. (pa|) reads: 



p(0=p i f(0 ; /(0 



and Eq. (|ll|), after integration, takes the explicit expression: 
e s 

^(0 



(58) 



(59) 



1+C 1+3 

where the truncated, scaled radius, defined by Eq. ([?]), is the analogon of the concentra- 
tion, c, defined in Navarro et al. (1997). 



The combination of Eqs. (03), (OH), (17), (63), and <M) yields 



M(£) = 12Af f 
p(0 = 12pt 1 



(60) 
(61) 



which represent the mass and mean density enclosed within a generic, isopycnic surface. 
It is worth noting that the profile parameter, Vp = p(S)/p' , is linked to the dimensionless 
density, 5 C , defined in Navarro et al. (1997), via the relation 5 C — 800/Vp. 



The profile parameters, v mas , v inri and v se i, due to Eqs. (14b), (15b), ( |21b| ), and (|59|), 
after integration, take the explicit expressions: 



^ / mo,.s — 12 



V se l = 36 



5 ' 3 V + 3 :- 6 » + e ln( i + 5 ) 

3(2 + 3) -2(1 + 3) ln(l + S) 



(62) 
(63) 
(64) 



(1 + 3)2 

the velocity profile related to centrifugal support along a radial direction, by use of 
Eqs. (H|) and ©, is: 

1/2 

; (65) 

and the profile parameter, v anm , due to Eq. (|58|), after integration, takes the explicit 
expression: 





'lln(l 


+0 




K)" 


1/2 


['7(01 




In 2 


-1/2 






[ ^ \ 



43 



S(2+ 3-21n ( l + S) 



i + a 



(66) 



which completes the specification of the physical parameters of interest, for the case 
under discussion. 

For further analysis of the properties of mass distributions with NFW density profiles 
(which is outside the aim of the current paper) see e.g., Lokas & Mamon (2001). 
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3.2. The MO A density profile 

The MOA density profile corresponds to the choice of exponents (a, j3, 7) = (3/2, 3, 3/2) 
in Eq. (|]). Accordingly, Eq. ( pa] ) reads: 

2 



p(0=P f /(0 ; f(0 



£3/2( 1 + £3/2) ' 

and Eq. (|i~l|), after integration, takes the explicit expression: 

4 (I icl/2^2 o O" 1 / 2 -! 

4 (i + £i/2)2 8 2 £i/2 _j 

the combination of Eqs. ©, (0), ©, and © yields: 



M(0 = 4Af t ln(l + £ 3/2 ) 
1 



p(0=V 7I ln(l + C 3/2 ) 



(67) 



(68) 



(69) 
(70) 



which represent the mass and mean density enclosed within a generic, isopycnic surface. 



The profile parameters, v mas and fi nr , due to Eqs. (141), (15b), (21b), and (|68|), after 
integration, take the explicit expressions: 



41n(l + ~ 3 / 2 ) ; 
13I/2 4. ^2 

4 



• _ 4 2S 1 / 2 -! 
■"inr - ' + ^ + -j= arctg -= 



-ln(l + S V V 2 ln(i _ sl/2 + s) + ^ 



(71) 



(72) 



the profile parameter, v se i, could also be analytically integrated^, but the result is ex- 
ceedingly long and, in practice, unserviceable. For this reason, it has to be calculated 
numerically. 

The velocity profile related to centrifugal support along a radial direction, by use of 
Eqs. (1253) and ©, is: 



nn(i + e 3/2 ) 



In 2 



,1/2 


\r)(0~ 


1/2 




[ . 





(73) 



and the profile parameter, v anm , due to Eq. (|67|) , after integration, takes the explicit 
expression: 



Van.m. — 2^ 



2 2S 1 / 2 - 1 2, , 1/9s 

arctg— — + - Hl + a W) 



(74) 



which completes the specification of the physical parameters of interest, for the case 
under discussion. 



2 By visiting the internet site: 



' http://INTEGRALS.WOLFRAM.COM/INDEX. CGI'' 
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The H density profile corresponds to the choice of exponents (a, 0, j) = (1, 4, 1) in Eq. (|l|). 
Accordingly, Eq. ( f^a]) reads: 

(75) 

(76) 

(77) 
(78) 



p(0=pV(0 ; f(0 = e(1 + 3 ; 

and Eq. (|TT|), after integration, takes the explicit expression: 

^(0 



(i + 2 (i + s) 2 ' 

the combination of Eqs. @, ©, (|l|), @, and © yields: 

e 



M(£) = 
#0 - 12p 



(i + 2 

fiJL 



£ 3 (i + 2 ' 

which represent the mass and mean density enclosed within a generic, isopycnic surface 



The profile parameters, v mas , v inr , and v se i, due to Eqs. (|l4b|) , ( |15b[ ), ( |21b| ), and (|76|), 
after integration, take the explicit expressions: 

"2 



12 



12 



^sei = 12 



(1+5)2 
' 5(2 + 5) 
1 + 5 
5 3 (4 + 5) 



21n(l + S) 



(79) 
(80) 
(81) 



(1 + 5)4 > 

the velocity profile related to centrifugal support along a radial direction, by use of 
Eqs. (HI) and ©, is: 



2£ 



1/2 



»?(0 



1/2 



(82) 



and the profile parameter, u anm , due to Eq. (|75|) , after integration, takes the explicit 
expression: 

5(2 + 35) 



45 



21n(l + 5) 



(83) 



(l + S) 2 

which completes the specification of the physical parameters of interest, for the case 
under discussion. 



3.4. Comparison with observations and simulations 

Self-similar, universal density profiles are characterized by two independent parameters, 
i.e. a scaling density, p\ and a scaling radius, r\ according to Eq. (Q). In some (e.g., 
NFW and MOA density profiles) but not all (e.g., H density profile) cases, the density 
must be null outside a truncation radius, to avoid an infinite, total mass. The truncation 
radius cannot exceed the tidal radius, due to the presence of neighbouring objects, and 
may safely be put equal to the virial radius (e.g., Cole & Lacey 1996; Navarro et al. 
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1997; Fukushige & Makino 2001; Klypin et al. 2001). In fact, the inner, denser regions of 
density perturbations collapse and virialize first, while the outer, less dense regions are 
still expanding or falling in. 

High-resolution, dark matter halo simulations from hierarchical clustering in either 
CDM or ACDM scenarios, are known to be consistent with both NFW and MOA density 
profiles (e.g., Fukushige & Makino 2001; Klypin et al. 2001). Then the particularization 
of homeoidally striated, density profiles with a central cusp discussed in Sect.||, to the 
special cases NFW and MOA, makes a useful tool for investigating some basic properties 
of large-scale celestial objects, such as galaxies or clusters of galaxies. 

With regard to elliptical galaxies, some properties of the stellar subsystem are repro- 
duced, to an acceptable extent, by H density profiles, in agreement with the observations 
(e.g., Hernquist 1990; Holley-Bockelmann et al. 2001). It is worth remembering that H 
density profiles closely approximate the de Vaucouleurs r 1 / 4 law for elliptical galaxies. In 
addition, recent observations (Geller et al. 1999; Rines et al. 2001) allow for the Coma 
cluster of galaxies both NFW and H density profiles, to represent the whole mass dis- 
tribution within about 10/z -1 Mpc, for a virial radius equal to 1.5/i _1 Mpc. The results 
of the current paper were applied to the Coma cluster of galaxies in a previous attempt 
(Caimmi 2002), and for this reason our attention shall be restricted to elliptical galaxies, 
where the boundaries of the dark matter halo and the baryonic ellipsoid may safely be 
idealized as similar and similarly placed. 

4. Application to elliptical galaxies 

According to current cosmological scenarios (e.g., Navarro et al. 1997), density pertur- 
bations at recombination epoch (z 1400) initially expand with the universe, turn 
around, collapse, and finally virialize (at least in their inner and denser regions). Just 
after virialization has been attained, baryonic and (dissipationless) non baryonic matter 
are expected to fill the same volume, and described by the same density profile. At late 
times the situation changes, as energy dissipation into heat within the gas subsystem 
makes it undergo further contraction, whereas the (non baryonic) dark matter does not 
appreciably change in extension. 

Virialized density perturbations, such as elliptical galaxies and clusters of galaxies, 
may safely be idealized as two homeoidally striated, similar and similarly placed, den- 
sity profiles with a central cusp. Our attention will be focused here on (giant) elliptical 
galaxies. 

4.1. General considerations and main assumptions 

A recent investigation performed on an optically complete sample of 42 (giant) elliptical 
galaxies, for which X-ray gas temperatures and central stellar velocity dispersions were 
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determined (Davis & White 1996), has shown evidence that (giant) elliptical galaxies 
contain substantial amounts of dark matter in general (Loewenstein & White 1999, here- 
after quoted as LW99). Accordingly, more than about 20% and 39%-85% of the total 
mass within one and six optical radii, respectively, is in form of (non baryonic) dark mat- 
ter, depending on the stellar density profile and observed value of X-ray gas temperature 
and central stellar velocity dispersion. The comparison between the velocity dispersion 
distributions for the dark matter and the stars, assuming isotropic orbits, shows that the 
dark matter is dynamically "hotter" than the stars, by a factor 1.4-2 (LW99). 

The above investigation cannot be considered as conclusive in favour of the existence 
of dark matter haloes hosting (giant) elliptical galaxies. In fact, it has been pointed out 
that the attenuation (in particular, the scattering) by dust grains has the same effect 
on the stellar kinematics as a dark matter halo (Baes & Dcjonghc 2001). On the other 
hand, current cosmological scenarios (CDM, ACDM) predict dark matter haloes hosting 
elliptical galaxies, as well as spiral galaxies, in the latter case supported by empirical 
evidence (e.g., flat rotation curves well outside optical radii). For this reason, we assume 
that (giant) elliptical galaxies are also embedded within dark matter haloes. 

An analysis on the evolution of the physical properties of cosmological baryons at 
low redshifts (z ~ 5) has recently been performed (Valageas et al. 2002), which (i) yields 
robust model-independent results that agree with numerical simulations; (ii) recovers the 
fraction of matter within different phases and the spatial clustering computed by numer- 
ical simulations; (iii) predicts a soft X-ray background due to the "warm" intergalactic 
medium component, that is consistent with observations. The related baryon fraction 
in the present universe is found to be 7% in hot gas, 24% in the warm intergalactic 
medium, 38% in the cool intergalactic medium, 9% within star-like objects and, as a still 
unobserved component, 22% of dark baryons associated with collapsed structures, with 
a relative uncertainty no larger than 30% on these numbers. Then the amount of still 
undetected baryons is about one fifth of the total, one fourth of the observed baryons 
(intergalactic medium, stellar components, and hot gas), and at least twice the stellar- like 
component. 

Given a typical (giant) elliptical galaxy, a natural question is to what extent the 
distribution of undetected baryons influences the "temperature" of both stellar and (non 
baryonic) dark matter. Towards this aim, we make the following main assumptions: 
(a) the stellar and the dark matter distributions are described by homeoidally striated, 
similar and similarly placed, H and NFW density profiles, respectively; (b) undetected 
baryons trace dark matter haloes; and (c) the virial theorem in tensor form holds for 
each subsystem. Let us discuss it briefly. 

Elliptical galaxies embedded within dark matter haloes, may safely be idealized as two 
homeoidally striated, similar and similarly placed, density profiles with a central cusp. 
According to recent investigations (e.g., LW99), viable representations for the outer, 
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non baryonic, and the inner, baryonic subsystem, are NFW and H density profiles, re- 
spectively. The above mentioned mass distributions were found to be self-consistent, in a 
parameter range of interest, with regard to the non negativity of the distribution function 
(LW99) by use of a theorem due to Ciotti & Pellegrini (1992). 

If undetected baryons in (giant) elliptical galaxies are present as hot gas, the gaseous 
subsystem is expected to be less concentrated than the stellar one, as in the Coma 
cluster of galaxies (e.g., Briel et al. 1992). In this respect, the simplest assumption is 
that undetected baryons trace the related, dark matter halo, i.e. are described by NFW 
density profiles. 

The typical velocity dispersion components, deduced by use of the virial theorem in 
tensor form, are global quantities, related to the potential-energy tensors of the subsys- 
tem as a whole, and so, by construction, independent of the specific orbital distribution of 
the particles. This important property, however, is also a weakness of the virial theorem 
in tensor form, when it is used to discuss velocity dispersion components measured in 
the central region of a galaxy. In fact, it is well known that the related values can be 
significantly different for structurally identical subsystems (and so characterized by iden- 
tical virial velocity dispersion components), due to different orbital structures (e.g., dc 
Zeeuw & Franx 1991). When using central velocity dispersion components, an approach 
based on Jeans equations (even though still questionable) is to be preferred (e.g., Ciotti 
& Lanzoni 1997; LW99). On the other hand, a comparison between the results obtained 
by use of either above mentioned methods, may provide additional support to both of 
them and/or useful indications on the nature of the problem under investigation. 

Strictly speaking, the central velocity dispersions (along the line of sight) in elliptical 
galaxies, which are deduced from observations, should be scaled to the corresponding 
typical values which make the virial theorem in tensor form hold (hereafter referred to, in 
general, as the virial velocity dispersions). Both observational evidence (e.g., Gerhard et 
al. 2001) and theoretical arguments (e.g., Nipoti et al. 2002) point towards the existence of 
dynamical homology in elliptical galaxies. In particular, a linear relation is found between 
a local parameter, averaged central velocity dispersion, and a global parameter, inferred 
maximum circular velocity, cto.i = (2/3)(v c ) ma;r (Gerhard et al. 2001). Accordingly, the 
central velocity dispersion components are expected to be proportional to the virial 
velocity dispersion components. Then it is assumed that the related proportionality factor 
is of the order of unity. 

In fact, typical peculiar velocity component distributions within (giant) elliptical 
galaxies show a maximum which is rapidly attained in the central region (at about 1 
kpc), and a slow decrease occurs moving outwards (no more than about 13% the maxi- 
mum at about 10 kpc), at least in the case of isotropic orbits; see e.g., LW99. Accordingly, 
both the central and the virial velocity dispersion components are expected to be of com- 
parable order, slightly less than the maximum of the peculiar velocity component distri- 
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bution. On the other hand, most elliptical galaxies are moderately radially anisotropic 
(e.g., Gerhard et al. 2001), and the related variation in central velocity dispersion (an 
increase for increasing <r^ and vice versa) is also expected to be moderate. 

4.2. Input parameters, specific assumptions, and results 

According to the above assumptions, a (giant) elliptical galaxy is idealized as two 
homeoidally striated, similar and similarly placed matter distributions, where the star 
and non baryonic subsystem are described by H and NFW density profiles, respectively. 
From this point on, the inner and the outer subsystem shall be labelled as *, D, instead 
of i, j, respectively. Let us suppose, at the moment, that undetected baryons are absent. 

Following LW99, we assign to the stellar subsystem a scaling radius r\ = (9/20)r e //, 
and a truncation radius i?* = 6 r e ff, where r e ff is the optical effective radius. Then the 
truncation, scaled radius, by use of Eq. (Q), is: 

R * 40 f OA \ 

-* = ^ = T ; (84) 

independent of the optical, effective radius. When the calculation of the scaling and the 
truncation radius is needed, the typical values: 

r eff = 5.04 ft^kpc ; (85) 
r\ = 2.268 ft- x kpc ; (86) 
fl* = 30.24 h^kpc ; (87) 

are used, where h is the dimensionless value of the present-day Hubble parameter, nor- 
malized to H = 100 km s _1 Mpc -1 . 

Turning our attention to the dark matter halo, and following again LW99, we assume 
for the truncation radius and the truncation, scaled radius, the typical values: 

R D = 256 h^kpc ; (88) 

Sz, = 10 ; (89) 
which, using Eq. (Q), yield the following, typical value for the scaling radius: 

r D = 25.6 fc-^kpc ; (90) 

in addition, the dark halo appears to be less concentrated than the stellar ellipsoid, i.e. 
Sj3 < 3* as expected from energy dissipation due to e.g., inelastic collisions between 
pre-stellar clumps. 

The mass ratio of non baryonic to baryonic matter related to a (giant elliptical) 
galaxy, equals the ratio of non baryonic to baryonic mean density within the volume of 
the system, p~d{Rd)/ P*(Rd), or the ratio of non baryonic to baryonic density parameter, 
flo/^b, as (e.g., Caimmi 2002): 



20 R. Caimmi and C. Marmo: density profiles with a central cusp 

where Qm = &d + &b is the (total) matter density parameter, and Qb = (0.05 =p 
0.01)(2/i) -2 fits to a good extent data on primordial nucleosynthesis (e.g., White & 
Fabian 1995). Accordingly, we assume the following value for the baryonic matter den- 
sity parameter: 

Q b = 0.0125 IT 2 ; (92) 

for a standard ACDM cosmological model with VIm = 0.3, Oa = 0.7, the combination of 
Eqs. © and © yields (Caimmi 2002): 

m D ,b = 24 h 2 - 1 ; (93) 

where 5 < mu,t < 23 as 0.5 < /i < 1. 

The mass of each subsystem, due to Eq. (^l|), keeping in mind that Mb = M* in the 
case under discussion, is: 

M. = j^— ; M D = ; (94) 

1 + rriD,* 1 + rriD,* 

where Mr = M* + Mb is the total mass, which has to be related to the virial mass, i.e. 

the mass within the virial radius of the corresponding density perturbation. 

Taking h = 2" 1 / 2 as fiducial value, and leaving n M = 0.3, n A = 0.7, Eqs. @, @, 

(H), ©, ©, and © yield: 

rj = 3.21 kpc ; R* = 42.77 kpc ; (95a) 
r D = 36.20 kpc ; R D = 362.04 kpc ; (95b) 
n b = 0.025 ; m D ,* = H ; (95c) 



and the combination of Eqs. (p5|), (|36|), (^5aj), and ( |95b| ) produces: 

y f = 11.29 ; y = 8.47 ; (96) 

the additional assumption that typical dark matter haloes hosting (giant) elliptical galax- 
ies have masses equal to 5 1O 12 M0, using Eqs. ([)i]) yields: 

M* = 4.54 lO n M ; M D = 5 10 12 M Q ; M T = 5.45 10 12 M Q ; (97) 

which completes the definition of the input parameters. 

The virial theorem in tensor form may be written for each subsystem separately (e.g., 
Brosche et al. 1983; Caimmi & Secco 1992) as: 

[(E U v)vir]pq + 2[(E u )kin]pq = ; u = *,D ; v = D, * ; (98) 

where [(E u )kin]pq represents the kinetic-energy tensor of uth subsystem. Let us define a 
typical velocity, v u = [(v u )i, (v u ) 2l (vu)a], as: 

[(E u )kin] P q = ^ vg M u (v u )p(v u ) q ; u = *,D ; (99) 
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then the combination of Eqs. fll4a| ), (|57|), (|98|), and (|99|) yields: 

1 K)p(<4)l , y. f v n n , 1m v 

: 1 = (Vu)mas(Vuv)vir ] U = *,D ; t> = L>, * ; (100) 

and the additional assumption of dominant random motions implies a velocity dispersion 
along the line of sight, (cr u ) p — (v u ) p in the special case where a principal axis points 
towards the observer, or using Eq. ( |100| ): 

1/2 

; u = *,D ; v = D,* ; (101) 



(cit)p 



| \yu)mas\yu 
(<lu)l 



(ep 2 £p3-4p) 1/2 

which has to be compared with its counterpart deduced from the observations. It is 
apparent that, in the absence of rotation, isotropy of peculiar velocity distribution occurs 
only for spherical configurations (E0), where £ P 2£p3^4p = 2/3, while the larger extent 
of anisotropy is attained for the most flattened (allowed) configurations (E7), where 
e 2 i = e 3 i = 0.3; Ai m 0.19; A 2 = A 3 « 0.90; oi > a 2 > a 3 . 

It is worth remembering that the axis ratios and the major semiaxis, which appear in 



Eqs. fllOOj ) and (101), are related to the intrinsic configuration. In doing the calculations, 
the scaling and truncation radius will be taken as representative of the major semiaxes 



of the corresponding, isopycnic surface. Using the results found in Subsects. 3.1 and 3.3, 



together with Eqs. fl95|), (|9q), (|97|), and ( 101 ), the velocity dispersion along the line of 
sight (assumed to coincide with the direction of a principal axis) takes the value: 

_ os7 i™ . ( fJl) )p _ oos i™ o-i 



h-^)* = 287 km 8 ; <MF~> = 225 km 8 ; <102) 

with regard to the stellar subsystem, it is consistent with a typical value er* = 250 km s _1 
(e.g., LW99; Treu et al. 2001) except in the limiting case of (intrinsic) E7 configurations 
with the major axis pointing towards the observer. 

The above results hold provided the baryonic subsystem is mainly in form of stars. 
Let us take into consideration a different scenario, where a less concentrated gaseous 
subsystem than the stellar one is also present, as in the Coma cluster of galaxies (e.g., 
Briel et al. 1992), and assume the same mass distribution as in the non baryonic mat- 
ter (Caimmi 2002). Accordingly, a (giant) elliptical galaxy may safely be conceived as 
formed by an inner subsystem made of stars and an outer subsystem made of gas and 
non baryonic matter. As the amount of baryonic and non baryonic matter have to remain 
unchanged, the inner and the outer subsystem are less and more massive, respectively, 
than in absence of undetected baryons. Again, it is assumed that the related mass dis- 
tributions correspond to a H and NFW density profile, respectively. Then we define the 
mass ratio: 

m g „ = ^ ; (103) 

where M* and M g are the mass of the stellar and gaseous subsystem, respectively. Bearing 
in mind that M*+M g = Mb, where Mb is the total mass in baryons, the following relations 
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are easily derived: 



M * = 1 , 1 ; M D + M g =[l+ ) M D ; (104a) 



Mb 

m D ,* + rri g ^ = (1 + m g ^)m D ,b + m g ^ ; (104b) 



and the repetition of the above procedure with M* = Mb, Mn, rrio,* — niD.b, replaced 
by M* = Mb — M g , Mo + M g , mo.* + mg,*, respectively, allows the specification of the 
input parameters. 

Some output parameters, such as profile factors involving both the inner and the 
outer matter distribution, and velocity dispersions along the line of sight, depend on the 
mass of each subsystem, and then on the mass ratio of gaseous to stellar subsystem, m g ^ 1 



via Eqs. (104)- The following results are found 
(o-*) p /(km s _1 ) 



(e P 2e P 3^p) 1/2 
(g D ) p /(km s' 1 ) 



= 246,223,208,198,190,184 ; (105) 
= 217, 213, 211, 210, 208, 208 ; (106) 



(e P 2e P 3^p) 1/2 

m g ,* = 0.5,1.0,1-5,2.0,2.5,3.0 ; (107) 

where m ffi * has been considered as a free parameter, in the range of interest. 

It is apparent that the outer subsystem (gas plus non baryonic matter) is dynamically 
"hotter" than the inner, provided the gas mass fraction exceeds the star mass fraction 
by a factor of about 3/2. The velocity dispersion along the line of sight decreases as m S) » 



increases, and a similar trend occurs, via Eqs. (101), (103), and (104), as the total mass 
decreases. 

On the other hand, the comparison (under the assumption of isotropic orbits) between 
the peculiar velocity component distributions related to the stellar ellipsoid and the dark 
matter halo, shows that the latter is dynamically "hotter" than the former (LW99). The 
reasons for this discrepancy may be due to (i) the different method used, with respect 
to the virial theorem in tensor form, in determining the "temperature" of the stellar 
ellipsoid and the dark matter halo (unfortunately, the derivation of the peculiar velocity 
component distribution is not outlined in LW99), and (ii) the different values of the mass 
ratio between the outer and inner subsystem in absence (mo,* = H) and in presence 



(ti_d.* = 29), according to Eqs. (104), near equally "hot" stellar ellipsoid and dark 

matter halo) of undetected baryons: values mo ,* + fng,* > 29 would allow dynamically 
"hotter" dark matter haloes, apart from the nature of their constituents. 

According to a recent investigation (Valageas et al. 2002), the ratio of undetected 
baryons associated with collapsed structures to star-like objects, in the present universe, 
attains a value which is close to 2. If undetected baryons and stars are present to a 
similar extent within (giant) elliptical galaxies, then m g .* ps 2 and dark haloes, according 



to Eqs. (105), (106), and (107), are dynamically "hotter" than stellar ellipsoids, in the 
case under discussion. In this view, both the observation that the temperature of the 
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extended hot gas exceeds the central stellar temperature, and the fact that the non 
baryonic matter is dynamically "hotter" than the stars (e.g., LW99), are a reflection of 
the presence of undetected baryons, which trace the dark halo and are about twice as 
massive as the stellar ellipsoid. 

Though (giant) elliptical galaxies host a surprisingly large amount of interstellar dust 
(up to a several ten million solar masses), most of it believed to be distributed diffusely 
over the galaxy (e.g., Baes & Dejonge 2001), this is still a negligible fraction of the total 
mass of the baryonic subsystem, such that a substantial amount of interstellar baryons 
within (giant) ellipticals, if it really exists, still has to be detected. 

5. Concluding remarks 

A general theory of homeoidally striated ellipsoids (e.g., Roberts 1962; Chandrasekhar 
1969; C93), where no divergence occurs in the density profile, has been adapted to cuspy 
density profiles. An explicit calculation of the related, physical parameters implies the 
specification of the density profile, which is equivalent to the knowledge of: (i) the func- 
tional dependence of a scaled density, / = p/p* , on a scaled radius, £ = r/r'; (ii) a 
boundary condition, i.e. /(I) = 1; (hi) two independent parameters, i.e. a scaling den- 
sity, p' , and a scaling radius, r' ; (iv) a truncated, scaled radius, 2. The latter requirement 
is due to the fact, that the systems under consideration exhibit a null density at an in- 
finite radius where, on the other hand, the total mass may attain a divergent value. In 
addition, an infinity of density profiles in the physical space, (Orp), is represented by a 
single density profile in the abstract space, (0£/), for any selected choice of exponents, 
(a, P, 7), appearing in Eq. @). 

Potential-energy tensors involving both one and two, homeoidally striated density 
profiles, where the boundaries are similar and similarly placed, have been expressed in 
terms of integrals on the mass distribution. Explicit calculations have been performed 
for both NFW and MO A density profiles, which satisfactorily fit the results of high- 
resolution simulations for dark matter haloes (e.g., Fukushige & Makino 2001; Klypin et 
al. 2001), and for H density profiles, which closely approximate the de Vaucouleurs r 1 / 4 
law for elliptical galaxies (e.g., Hernquist 1990; Holley-Bockelmann et al. 2001). 

The virial theorem in tensor form, related to a two-component system, has been 
expressed for each subsystem, and applied to giant elliptical galaxies. The predicted 
velocity dispersion along the line of sight, in the limiting case where a principal axis 
points towards the observer, has been found to be consistent with the data except for 
(intrinsic) E7 configurations, when the major axis points towards the observer. 

The suggestion that dark matter haloes host an amount of undetected baryons as 
massive as about twice the stellar subsystem (Valageas et al. 2002), together with the 
assumption that undetected baryons trace non baryonic matter therein, has produced two 
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main consequences, namely (i) predicted velocity dispersions along the line of sight are 
lower than in absence of undetected baryons, and (ii) dark matter haloes are dynamically 
"hotter" than stellar ellipsoids, the transition occurring when the amount of undetected 
baryons is about one and a half times the stellar subsystem. 

In this view, both the observation that the temperature of the extended hot gas 
exceeds the central stellar temperature, and the fact that the non baryonic matter is 
dynamically "hotter" than the stars (e.g., LW99), are a reflection of the presence of 
undetected baryons, which trace the dark halo and are about twice as massive as the 
stellar ellipsoid. 
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Appendix A: A procedure for fitting simulated density profiles 

Simulated, spherically averaged, virialized dark matter haloes, allow the knowledge of 
the density profile, the virialized mass, and the virialized radius, where the last two may 
safely be taken as representative of the whole mass and the whole radius, respectively 
(e.g., Cole & Lacey 1996; Navarro et al. 1997; Fukushige & Makino 2001; Klypin et al. 
2001). A plot of the density profile on the logarithmic plane, (0 log £ log /) , necessarily 
implies use of dimensionless coordinates, defined by Eqs. ( pa|) and (^), namely: 

/co = 4 ; £=-7 ; ( A -v 

p> r< 

where the choice of the scaling density, p\ and the scaling radius, , is, in principle, 
arbitrary. A different choice of the above mentioned parameters, say from (jy,r*) to 
[(//) , ( r ') ]> makes simulated points shift by a factor, log[(/?') ' / p'], along the log / axis, 
and by a factor, log[(r'y /r'], along the log£ axis. 

In fitting simulated density profiles, our attention shall be restricted to the family of 
functions (e.g., Hernquist 1990; Zhao 1996): 
A [3 — 7 

/(e) = r(i+e)* 5 x = ^r ; (A - 2) 

which has been deduced from Eqs. ([!]) and (A.l), provided the boundary condition: 
p(r+)=pt/(i) = p t ; /(1) = i ; J=p{l) ■ (A.3) 
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holds, according to Eq. (J2a|). The combination of Eqs. (|A.2|) and (A. 3) yields 



f(l)=A2-* = l ; A = 2* ; (A.4) 



and Eq. ( |A.2| ) takes the final form: 
which is equivalent to: 

log/ = xlog2- 7 log^-xlog(l + D ; (A.6) 



due to Eq. (A.3), the three-parameter curve represented by Eq. (A.6) remains unchanged 



in the logarithmic plane (0 log £ log/), for different choices of the scaling parameters, 



(p',r*). On the other hand, due to Eq. (A.l), any point (£, /) on the curve represents 
different points (£r' t fp') on (one or more) physical density profiles, for different choices 
of the scaling parameters, (jy,r'). 

In the limit of negligible values of the independent variable, £, with respect to unity, 



Eq. (AJ3) reduces to: 



log/ = xlog2-7log£ ; e«l ; (A.7) 

which represents, in the logarithmic plane, a straight line with slope equal to —7 and 
intercept equal to xlog2. 

In the limit of preponderant values of the independent variable, £, with respect to 



unity, Eq. (AJ3) reduces to 



log/ = X log2-/?log£ ; e»I ; (A.8) 

which represents, in the logarithmic plane, a straight line with slope equal to —j3 = 
— (7 + \a) a n d intercept equal to xlog2. 

The straight lines under discussion have coinciding intercepts i.e. they meet on the ver- 



tical axis and, in addition, represent the asymptotes of the curve, expressed by Eq. (A.6) 



The special cases related to NFW and MOA density profiles, are plotted in Fig. A.f . 

The above results hold for a > 0. The case a < makes the asymptotes change one 
into the other. The limiting case a = makes the asymptotes coincide i.e. the curve 
reduces to a straight line. 

A change of variables: 

z = log£ ; y = logf ; (A.9) 
translates Eq. flA.fiD into the equivalent form: 



V = X log 2 - -fx - x log[l + exp 10 (ax)] 



(A.10) 
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Fig. A.l. Representation, in the logarithmic plane, of NFW (left) and MOA (right) 
density profiles (full curves), with their asymptotes (dashed lines). In both cases, — 1 < 
log£ < +1 and -3 < log / < +3. 



where, in general, exp u (x) — u x , and exp(x) = e x , according to the standard notation. 
The derivatives up to the third order are: 



-0- 



Xa 



dy 
dx 

— ^ = -\a ln(10) 
4jj = -xa 3 ln 2 (10)exp 10 (a:a;) 



1 + exp 10 (aa;) 

exp 10 (aa;) 



exp 10 (aa;)] 2 

1 — exp 10 (aa;) 



[1 +exp 10 (aa;)] 3 



and the particularization of Eq. (A. 11a) to x — i.e. ^ = 1 or r = r>, reads 



(A.lla) 
(A.llb) 
(A.llc) 



(A.12) 

in the special case of NFW density profile, the right-hand side member reduces to -2, as 



pointed out by Bullock et al. (2001). On the other side, Eq. (A.llc) shows that the third 
derivative equals zero at the same point, x = i.e. £ = 1 or r = r'. Accordingly, density 
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profiles expressed by Eqs. (A.l) and (A. 2), are characterized by a maximum variation in 
slope, in the logarithmic plane (0 log£ log/), at that point where x = i.e. £ = 1 or 



and the slope is expressed by Eq. ( A . 1 2| ) 



Bearing in mind that fixed points on a physical density profile change their position on 
the scaled density profile, for different choices of the scaling parameters, {p\ r'), according 



to Eq. (A.l), a possible procedure in fitting density profiles of simulated, virialized, dark 
matter haloes, could be the following. 



(i) Assume that self-similar, universal density profiles, expressed by Eq. flA.q ), are rep- 
resentative of virialized, dark matter haloes (e.g., Cole & Lacey 1996; Navarro et al. 
1997; Fukushige & Makino 2001; Klypin et al. 2001). 

(ii) Select two special choices of scaling parameters, [(pt)~, and [(p^) + , ) + ], 
which make simulated points well shifted on the left and on the right, respectively, 
i.e. (expected to be) sufficiently close to either asymptote. 

(iii) Select the best pair of deduced asymptotes, which meet together on the vertical 
axis. 

(iv) Calculate the slope and the intercept of either asymptote. 



(v) Calculate the exponents, (a,(3,y), via Eqs. (A. 2), (A. 7) and (A 



In substance, the above procedure relies on the possibility, that a least-square fit 



related to a curve, represented by Eq. (|A.6| ), may be reduced to a least-square fit related 



to a straight line, represented by either Eq. ( A. 7 ) or Eq. (A 



We hope that, in fitting self-similar, universal density profiles of simulated, virialized, 



dark matter haloes, with mass distribution expressed by Eq. (A.l), the procedure will be 
explained in detail. 



Appendix B: Connection between NFW and MOA density profiles in fitting 
simulated haloes 

A generic density profile, belonging to the family expressed by Eq. (|l|), is defined by 
three exponents, (a,/3, 7), and is characterized by two independent parameters e.g., a 
scaling density, p\ and a scaling radius, r', or any other pair of independent parameters. 
In particular, the total mass within the virial radius, M = M20O) an d a dimensionless 
parameter, S, may be used, and their expression as a function of and r' are derived 
by fitting numerical simulations (Fukushige & Makino 2001, hereafter quoted as FM01). 

Then the problem is how to compare simulated density profiles with different, theoret- 
ical ones of the kind under consideration. In doing this, our attention shall be restricted to 
(a) high-resolution simulations from FM01, and (b) NFW and MOA theoretical, density 
profiles, but the method is quite general. 
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With regard to NFW density profiles, the best fit to numerical simulations prescribed 
by FM01: 

m= (£/2)[l + (£/2)P ; (R1) 
is inconsistent with the curve plotted in their Fig. 18 (hereafter quoted as Fig. FM18), 
probably due to some printing errors. In fact, the comparison with their Fig. 17 (hereafter 



quoted as Fig. FM17), which shows the counterpart of Eq. (B.l) related to MOA density 
profile: 

/(£) = £3/2( 1+ £3/2) i ( R2 ) 

allows the following conclusions: (i) simulated density profiles remain unchanged, when 
both the above mentioned plots have same origin and same scale; and (ii) both NFW 
and MOA density profiles, represented therein, attain the same value as £ — * +oo. 



With these constraints, Eq. (BA) takes the form: 

m = (TOf ; (R3) 

where k\ is a parameter which can be calculated from the knowledge of the coordinates 
of an arbitrary selected point on the NFW density profile, plotted in Fig. FM18. 



Following FM01, NWF and MOA density profiles are normalized as in Eqs.(B.l) 



and (B.2), respectively, and plotted in Figs. FM17 and FM18, respectively. On the other 



hand, Eq. (B.l) appears to be inconsistent with the curve plotted in Fig.FM18. To get a 
deeper insight into this problem, let us start with the obvious relations: 

p** = — 7 — ; (b.4) 

\P J MOA 

; (B.5) 



p 


P 


(p')moa 


„ (p')moa 


(p')nfw 


(p')moa 


(p')nfw 


\P JNFW 


r 


r 


( r ')MOA 


_ „ ( r ')\lOA 


( r ')NFW 


( r ')MOA 


( r ')iVPW 


( r ) NFW 



V JMOA 

where defined as in FM01, and (p')\ (r') f , have the same values as prescribed 

therein. 

Bearing in mind that any point, P(r**, /?**), maintains its coordinates passing from 
Fig. FM17 to Fig. FM18 and vice versa, the NFW density profile plotted in Fig. FM18, 
reads: 

k 1 k 2 



(&ir**)[l + (fcir**)] 5 

]„ - (T%£QA . /- _ 1 (p')nfv . <R(,i>< 
1 ~ (r') f ' h (P') f ' { ' 

V 1NFW 1 \P IMOA 

further inspection of Figs. FM17 and FM18 discloses the validity of the condition: 

-3/2/-. 3/2s_j 

lim r ** 1 + r ** > = 1 ; (B.7) 

r«^+oo k 2 r^ (1 + fcir»»)~ 2 

which, performing the related calculations, is equivalent to: 



(B.6a) 
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and the combination of Eqs. (B.6a) and (B.8) yields Eg. (JB . 3|) . The comparison of 



Eqs. flB.l| ), ( [B.2D , with (|5J), (|67j), respectively, yields: 

2(p f )MOA = {p')moa > (^)moa = (r'^uoA ; ( R9a ) 
Mp^nfw = kl{p')uo A ; ( rt )jv-FW = fc 1 - 1 (r')Mo^i ; ( R9b ) 
which translate the prescriptions from FM01, their Eqs. (13) and (14), into the following: 
, _ - - / M V' lg°Mg 



(Cp)atfh' = -^r ; (C p ) A fOA = 2Q ! (B.lOb) 

with regard to the scaling density, , and: 
/ /if \ 2 / 3 

rt = ar 1 ' 3 ^— J kpc ; (B.lla) 

(C r )jvFW = 2 lCT 2 / 3 ^" 1 ; (a)MOA = 2 1(T 2/3 ; (B.llb) 

with regard to the scaling radius, r\ where Mio = 1O 1O M0 and the values of {p*)nfw, 
(p )moa, are one fourth, one half, respectively, the value derived in FM01, due to the 
different normalization of the NFW, MOA, density profiles, adopted in the current paper. 
The value of the dimensionless parameter, 5, is considered to reflect an amplitude of the 
density fluctuation at turnaround and, for this reason, it can be thought of as constant 
during the evolution of a halo (e.g., Cole & Lacey 1996; Navarro et al. 1997; FM01). 
From the standpoint of top-hat, spherical density perturbation, it is related to both the 
mass and the peak height. 

At this stage, what still remains to be done is the specification of the NFW density 
profile related to the fit of simulated, dark matter haloes from FM01 plotted in Fig. FM18. 
Towards this aim, let P(f**,/5**) be a selected point on the curve plotted in Fig. FM18. 
Accordingly, Eq. ( |B.3| ) reads: 

" p ** = Mi + (W)] 2 ; (R12) 

and the related, second-degree equation in k\ has the solutions: 

kf = r- X {[l ± (P**^*) 172 ]- 1 - 1} ; (B.13) 

where the negative solution has no physical meaning. Then the acceptable solution to 
Eq. fllTT^ ) is: 

k 1 =r^ 1 {[l-(p**rl) 1 / 2 )- 1 -l} ; (B.14) 

where it is intended that the square root is non negative. 

The results related to a pair of selected points on the curve plotted in Fig. FM18, are 



listed in Table B.l. The counterpart of Eq. (B.3) in FM01, probably due to some printing 
errors, is inconsistent with the results listed in Table [B.l| , as it would imply kf = 10, 



fci = 0.5, according to Eq. (B.l). Consistent, even if non acceptable, choices would be 
either kf = 8, ki = 2, or kf = 10, ki = 10 1/3 = 2.1544. 
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f>„ 




ki 


7 2 

kl 


kf 


1000 


0.005 


2.26 


5.11 


11.56 


500 


0.01 


2.29 


5.23 


11.91 



Table B.l. Specification via the parameter, k\, of the NFW density profile related to the 
fit of simulated, dark matter haloes from FM01, plotted in Fig. FM18. The coordinates 
of the two points, P(f**,p**), have been deduced from Fig. FM18. 



